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ABSTRACT 

We present the discovery of a substellar companion on a wide orbit around the —2.5 Mq star HIP 78530, 
which is a member of the 5 Myr-old Upper Scorpius association. We have obtained follow-up imaging over 
two years and show that the companion and primary share common proper motion. We have also obtained 
JHK spectroscopy of the companion and confirm its low surface gravity, in accordance with the young age 
of the system. A comparison with Drift-Phoenix synthetic spectra indicates an effective temperature of 
2800 ± 200 K and a comparison with template spectra of young and old dwarfs indicates a spectral type of 
M8±l. The mass of the companion is estimated to be 19-26 Mj up based on its bolometric luminosity and the 
predictions of evolutionary models. The angular separation of the companion is 4.5", which at the distance of 
the primary star, 156.7 pc, corresponds to a projected separation of —710 AU. This companion features one of 
the lowest mass ratios (—0.009) of any known companion at separations greater than 100 AU. 

Subject headings: stars: pre-main sequence — stars: low-mass, brown dwarfs — stars: formation — planetary 
systems 



1. INTRODUCTION 

The recent direct imaging discoveries of very low mass 
substellar companions to stars - especia lly those orbiting the 
young stars 1RX S J160929. 1-21052 4 (lLafreniere et al.ll2008l 



2010), HR 8799 (Marois et al. 2008), Fomalhaut (Kalas et al. 
2008), AB Pic (IChauvin et alj|2005j>. HN Peg dLuhman et all 
20071) . and CT Cha (ISchmidt et alj|2008b - have raised new 
questions about the formation mechanisms of planets and 
brown dwarfs. These companions have separations of tens 
to several hundred AU and mass ratios <0.02 relative to their 
primaries. Several theoretical ideas are under discussion to 
explain such properties, but they all run into some difficul- 
ties. Core accretion models favor formation of giant plan- 
ets cl ose to the 'snow line', i.e., at <10 AU (Pollac k et al.l 
1996). Gravitational instability could produce massive plan- 
ets at large radii, but would r equire unusually large an d mas- 
sive circumstellar disks (e.g. Vorobvov & Basu 2010). Frag- 
mentation of pre-stellar cores during colla pse does not lea d 
to such extreme binary systems easily (e.g. lBate et al.ll2003h . 
Thus, some theorists have recently proposed that dynamical 
instabilities within planetary systems that originally formed 
multiple giant planets could scat t er some of them t o large sep- 
arations (Scharf & Menou 2009t lVeras et al-O OOSf). However, 
the more massive planets are less likely to reach the widest 
orbits. 

Here we report on the direct imaging discovery, common 
proper motion confirmation and multi-band spectroscopy of a 
-23 Mj up companion seen -700 AU from HIP 78530, which 
is a —2.5 Mq star (spectral type B9V) in the Upper Scorpius 
young association. The separation of this new companion 
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places it among the widest known substellar companions to 
stars and its extreme mass ratio - among the lowest currently 
known - is comparable to those of directly imaged planets, 
even though its mass is well above the deuterium-burning 
threshold. Together with the companions mentioned above, 
this new substellar companion presents a good challenge to 
all formation scenarios and contributes to blurring the distinc- 
tion between giant planets and brown dwarfs even further. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Imaging 

The discovery presented in this paper was made as part of a 
direct imaging search for new stellar and substellar compan- 
ions around about 90 stars in the Upper Scorpius region. The 
overall target sample was built b y randomly selecting, f rom 
the list of Upper Scorpius stars in ICarpenter et al.1 (12006ft . an 
equal number of stars in each of five equal logarithmic mass 
bins over the range —0.15-5 M Q ; the spectral types range 
from B O to M5. The observ ations were made using the NIRI 
camera (Hodarjp et al. 2003) and the ALTAIR adaptive optics 
system dHerriot et al.l 120001) at the Gemini North Telescope. 
Except for a few faint targets, the target stars themselves were 
used for wave front sensing. The field lens of ALTAIR was 
used to reduce the effect of anisoplanatism. The first epoch 
imaging of HIP 78530 was done on 2008 May 24 in the nar- 
row band filter ^c 0nt i nuurn centered at 2.0975 /im. For sky sub- 
traction we used 5 dither positions corresponding to the cor- 
ner and center of a square of side 10". At each position we 
obtained one co-addition of twelve 0.5 s integrations in fast, 
high read-noise mode, followed by one single 10 s integration 
in slow, low read-noise mode. At each position this provides 
an unsaturated image of the target star and a much deeper 
image of the field that can be readily spatially registered and 
scaled in flux. The full-width-at-half-maximum for this data 
set is 0.075" and the Strehl ratio is 0.34. 

The initial image of HIP 78530 revealed an interesting faint 
nearby source which, based on a c o mpari s on wi th observa- 
tions made by iKouwenhoven et ail d2005l [2007 ), was very 
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likely to be a true co-moving companion. To confirm the com- 
mon proper motion of this candidate and get additional pho- 
tometry measurements, follow-up imaging observations in J, 
H and K 1 were obtained on 2009 July 2 using the same in- 
strument as before. To avoid saturation in these broad fil- 
ters, we used a512x512 subarray, and kept the strategy of 
a short, unsaturated image immediately followed by a longer, 
saturated image at each of 5 dither positions. For each wave- 
length, three sky frames were obtained at offsets of > 15". 
The short exposure times were 0.055 s with 40 co-additions, 
while the long exposures were 6 s, 5 s, and 4 s with one co- 
addition in J, H and K', respectively. To increase the sig- 
nificance of our common proper motion confirmation further, 
another follow-up observation was made on 2010 August 30 
using the filter, full frame, and one co-addition of twelve 
0.6 s integrations followed by a single 10 s integration at each 
of five dither positions. Other follow-up data were acquired 
in spring and summer 2010 but owing to an overseen differ- 
ence in the instrumental setup, these data suffer from large 
systematic astrometric errors and are not used here. 

The imaging data were reduced using custom IDL routines. 
For the 2008 and 2010 imaging data, a sky frame was con- 
structed by taking the median of the images at all dither po- 
sitions after masking out the regions dominated by the tar- 
get's signal, while for the 2009 imaging data, the median of 
the three sky frames was obtained. After subtraction of this 
sky frame, the images were divided by a normalized flat-field. 
Then isolated bad pixels were replaced by the interpolated 
value of a third-order polynomial surface fit to the good pixels 
in a 7 x 7 pixels box while clustered bad pixels were simply 
masked out. Next the images were distortion corrected using 
the distortion solution provided by the Gemini staff0 Finally, 
the long-exposure images were properly scaled in intensity 
and merged with their corresponding short-exposure images. 
The color composite JHK image from the 2009 follow-up ob- 
servations is shown in figure Q] 

2.2. Spectroscopy 

The likely companion detected in the imaging data was 
very faint and, assuming it was a true companion, its mag- 
nitude suggested a substellar mass. Thus, to verify its late- 
type nature, we obtained spectroscopic observations of it in 
H, K and J on 2009 July 2, 2009 July 3, and 2009 August 
8, respectively, using the integral field spectrograph NIFS 
dMcGregor et al.ll2003l) with the ALTAIR adaptive optics sys- 
tem at the Gemini North telescope. The J grating was used 
with the ZJ blocking filter, the H grating with the JH block- 
ing filter and the K grating with the HK blocking filter. The 
spectral resolving power is ^6000 in J and ^5300 in both H 
and K. Given the large angular separation of the companion 
(~4.5"), the primary star is not visible in the 3" x 3" NIFS 
field-of-view. In each band, the companion was positioned 
near the center of the field-of-view and was dithered by 0.7" 
along a line between each of five exposures to enable good 
sky subtraction. The individual exposure times were 480 s, 
300 s and 240 s in J, H and K, respectively, in low read noise 
mode. After the J and H sequences, the A0V star HIP 79229 
was observed at a similar airmass for telluric and instrumen- 
tal transmission correction, while the A0V star HIP 73820 
was observed just before the K band sequence for the same 

5 The distortion correction used is centro-symmetric and is given by r = 
r' + 132 X 10~ 5 r' 2 , where r and r' are respectively the distortion-corrected 
and distorted radial pixel distances from the array center. 



TABLE 1 
Astrometric measurements 



Epoch 


Band 


P(") 


P.A. (deg) 


2008.3940 
2009.4998 
2010.6605 


"-cont 

JHK' 

r 2.QQ 
"cont 


4.529 ±0.006 
4.533 ±0.006 
4.536 ±0.006 


140.32 ±0.10 
140.35 ±0.10 
140.27 ±0.10 



purpose. 

The reduction of the NIFS data, up to the reconstruction of 
the data cube, was made using the Gemini IRAF pipeline. The 
steps covered in this pipeline are sky background subtraction, 
flat-field and bad pixel correction, spatial and spectral calibra- 
tion, and data cube reconstruction. The spatial sampling of 
the reconstructed data cube is 0.043" pixel" 1 . The instrumen- 
tal/telluric transmission correction and spectrum extraction, 
which could be done using the Gemini pipeline, were instead 
done using custom IDL. 

The center of each PSF was first registered to a common 
position in all spectral slices and all cubes of the sequence; 
the center positions were calculated by fitting a 2D Gaussian 
function. The spectrum of the source was extracted by sum- 
ming the flux in a circular aperture of diameter 5 pixels. The 
spectrum of the telluric standard was corrected for its spectral 
slope using a 9520 K blackbody curve, and any hydrogen line 
absorption was removed by dividing out a Voigt profile fit. 
The companion spectrum was divided by the standard spec- 
trum to correct for telluric and instrumental transmission. Fi- 
nally, the median of the 5 spectra was obtained. 

3. ANALYSIS AND RESULTS 

A summary of the properties of HIP 78530A and its com- 
panion, taken from the literature or derived in this section, is 
presented in Table [2] 

3.1. Common proper motion 

The relative position of the companion and primary was 
determined by fitting a 2D Gaussian function to each com- 
ponent. The pixel positions were converted to arcseconds us- 
ing the 2 1 .4 mas pixel scale of NIRI with ALTAIR and the 
field lens as indicated on the instrument web pag^5 the po- 
sition angle of the image was taken from the FITS header. 
The astrometric errors, 6 mas in separation and 0.1° in po- 
sition angle, were estimated in lLafreniere et alJ ((2010) using 
the position of background stars for a similar sequence of ob- 
servations; they are mostly dominated by residual distortion 
errors. We also note that the pixel scale of ALTAIR/NIRI 
with the field lens has not been extensively calibrated and 
that our values may be systematically different from measure- 
ments made using other instruments or telescopes, but they 
should be internallyconsistent. The astrometry measurements 
are given in Table [Hand shown in Fig. [2] in comparison with 
the changes expected over time for a distant background star, 
based on th e proper motion an d distance of the primary star 
(taken from Ivan Leeuwen 2007). Our measurements are con- 
sistent with common proper motion but inconsistent, at a level 
of ^6c, with the changes expected for a distant, motionless 
background star. This clearly indicates that the companion is 
co-moving with the primary. 

The accuracy of our astrometry measurements can be veri- 
fied using a faint background star detected both in 2008 and 

6 http : //www . gemini . edu/sciops /instrument s /alt air / field- lens- 
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2010. This source, ~ 10. 1 mag fainter than the primary, had a 
separation of 9.451" and position angle of 43.49° at epoch 
2008.3940, and 9.493" and 43.48°, respectively, at epoch 
2010.6605. These values, showing a ^7er change in sepa- 
ration and ~ lex in position angle, are fully consistent with 
expectations for a distant, motionless background star. 

The new companion we report in this paper was 
detected previously by jKouwenhoven et ail (120051) and 
iKouwenhoven et all(l2007l) in data obtained in 2000 and 2001 
with ADONIS at the ESO 3.6 m telescope. They report a sep- 
aration of 4.54" ±0.01" and position angle of 139.7° ±0.3°. 
They note that they were not able to conclude whether this ob- 
ject was a bound companion or a background object based on 
HK photometry, and they have not followed up this source 
further. Considering the uncertainties involved in compar- 
ing our astrometry measurements with those made using 
other instruments, as noted above, the astrometry reported by 
Kouwenhoven et al. is in very good agreement with ours and 
further supports our conclusion that this source is co-moving 



with the primary star. Indeed, were the source a background 
star, its separation should have changed by -0.1" and its po- 
sition angle by -3.3° between 2000 and 2010, while the mea- 
sured values show changes of +0.01" and +0.6°, respectively. 

3.2. Properties of the companion 

In this section, in addition to comparing our spectrum and 
photometry measurements with those of other objects, we also 
rely on comparisons with synthetic spectra to estimate the 
physical properties of the companion. For the comparison 
with mode ls, we adopted t he Drift-Phoenix atmosphere 
mode ls dWitte et all 120091; iHelling et al.1 120081; iDehn et alj 
120071) . which combine a detailed kinetic mo del of dust cloud 
formation with a radia tive transfer code (lHauschildt et alJ 
Il999t iBaron et al.ll2003l) . In contrast to all other atmosphere 
models which assume phase equilibrium between gas and 
cloud particles, Drift-Phoenix describes consistently the 
formation of a stationary cloud by homogeneous nucleation 
and grain growth/evaporation, including gravitational settling, 
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FIG. 2. — Measured (data points) separations (top) and position angles (bot- 
tom) between the primary star and its companion as a function of time. The 
solid lines show the expected position of a distant stationary background 
source over time, as calculated from the proper motion and distance of the 
primary star. 



TABLE 2 
Properties of HIP 78530 AB 



Value 



Parameter 


Primary Companion 


fXa cos 8 (mas yr ) a 


-11.38±0.59 


fM S (mas yr~') a 


-24.66 ±0.37 


Distance (pc) a 


156.7± 13.0 


Angular separation (") 


4.533 ±0.006 


Position angle (deg) 


140.3 ±0.1 


AJ (mag) 


8.13±0.05 


AH (mag) 
AK' (mag) 


7.44 ±0.03 


7.27 ±0.03 


AK™> (mag) 


7.26 ±0.03 


J (mag) b 


6.928±0.021 15.06±0.05 


H (mag) b 


6.946 ±0.029 14.39 ±0.04 


K s (mag) b 


6.903±0.020 I4.17±0.04 


J-K s (mag) 


0.025 ±0.03 0.89 ±0.06 


H-Ks (mag) 


0.043 ±0.04 0.22 ±0.06 


Spectral type 


B9V M8±l 


T at (K) 


~10500 c 2800±200 


log (L/Lq) 


-2.55±0.13 


Mass (Mq) 


~2.5 d 0.022 ±0.004 c 


Projected separation (AU) 


710 ±60 



" From van Leeuwen (2007). 

b From the 2MASS PSC and our contrast measurements. Differences in filter 

bandpasses may add up to 1-2% uncertainty. 

c From the spectral type based on the temperature scale of Sh erry et"aT] 



d From the models of D'Antona & Mazzitelli 1 1997). 

c From the models of Baraffe et al. 1 1998, 2002) and Burrow s~et"aT1 (19971) . 

element depletion, and convective element replenish ment 
dWoitke & HelUngll2Q0l l200l iHelling & Woitkell2006h . All 
models use d assume solar meta llicity, which is reasonable for 
Upper Sco dMohantv et al.ll2004l) . 

3.2.1. Temperature, surface gravity and spectral type 

The spectrum of HIP 78530B, binned down to a resolving 
power of 400, is shown in figure [3] in comparison with syn- 
thetic spectra from the Drift-Phoenix atmosphere models 
for various temperatures and surface gravities. As visible on 
the figure, for all effective temperatures shown, the spectrum 
of the companion is in better agreement with the lower sur- 
face gravity spectrum, as expected for a young object in Up- 
per Sco. This is particularly striking in the K band where the 



broad H2O absorption bands and the CO absorption band head 
depths are perfectly reproduced by the low gravity spectrum 
at r e ff = 2600 K. The low surface gravity is also apparent in the 
depth of the Kl lines in the J band. The effective temperature 
providing the best fit depends on the bandpass considered. As 
mentioned previously, a temperature of 2600 K provides an 
excellent fit in the K band. A temperature of 2600 K also 
provides a slightly better fit in the J band, both for the con- 
tinuum around 1.2 /im and the depth of the Kl lines. In the 
H band, however, the best fit is obtained for temperatures of 
^2800-3000 K. Taking these considerations into account, we 
adopt an effective temperature of 2800 ± 200 K. We have also 
done a chi-square minimization using the same set of model 
spectra, rather than a simple visual inspection, and reached 
the same conclusion. In particular, when the minimization is 
done in all bands simultaneously, the best fit is found for a 
temperature of 2800 K. 

Figure |4] shows the companion spectrum, at a higher reso- 
lution this time, in comparison with the spectrum of a known 
M8 dwarf member of Upper Sco as well as with a model 
spectrum with our effective temperature estimate of 2800 K. 
Among the Upper Sco brown d warfs identified and spectrally 
classified by lLodieue t al. (2008), the best fit to our spectrum 
is obtained for a spectral type of M8, which corresponds to an 
effective temperature of ^2700-2800 K, consistent with our 
previous estimate. Comparing our spect rum to those of field 
dwarfs from the IRTF spectral library ( R avner et al.l 120091) 
also results in a best fit for spectral types of M7-M9, although 
there are differences between the spectra owing to the lower 
gravity of the companion as noted earlier. The better agree- 
ment of our spectrum with those of other objects in Upper 
Sco, as opposed to field objects, provides further evidence that 
the new companion also belongs to the association. 

The companion spectrum is not only in good agreement 
with those of models and other Upper Sco objects within the 
individual bands, but it also shows a reasonable agreement 
in colors between the bands. Over a temperature range of 
2600-3000 K, the Drift-Phoenix models yield synthetic 
colors of 7-/7=0.50-0.55 and H -£=0.27-0.38, to be com- 
pared with the companion colors of J-H=0.67 ± 0.06 and 
H-K=0.22 ± 0.06. Assuming a s mall possible extinctio n 
of Ay = 0.5 mag toward this system (Carpe nter et al.l 2009). 
the extinction-corrected colors would be J - H^0.62 and 
U- .1 9. The colors of M8 dwarfs in Upper Sco from the 
sample of lLodieu et ail §008) are 0.58-0.65 and 0.38-0.57, 
respectively. Here, the agreement in J-H is good but our 
companion is significantly bluer in H - K. We note that the 
differences between the NIRI (used for the companion con- 
trast) and 2MASS (used for the primary star magnitude) filter 
bandpasses add an uncertainty of at most 2-3% on the colors 
of the companions, calculated from synthetic spectra. 

A good indicator of surface gravity is the equivalent width 
(EW) of the Kl atomic lin e doublets at 1.168/1. 177 /im 
and 1.243/1.253 /im (e.g. iMcGovern eTall 12004 . For 
HIP 78530B, the EWs of these four lines are, respectively, 
2.0 ±0.3 A, 3.0 ±0.3 A, 2.0 ±0.2 A, and 2.1 ±0.2 AQ In 
the Drift-Phoenix models, there is a clear monotonic trend 
that the EW in these lines increases with increasing surface 
gravity for a constant temperature, see figure[5] We have eval- 
uated the EW for models of 2800 K with surface gravity rang- 

7 To calculate the EWs, we used an integration interval of 10 nm centered 
on each line and fitted the continuum as a straight line using 40 nm-wide 
intervals on both sides of the line. 
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FIG. 3. — NIFS spectrum of HIP 78530B (black) compared with synthetic spectra from the DRIFT-PHOENIX atmosphere models (red and blue) of various 
effective temperatures and surface gravities. From the top to the bottom rows, the effective temperatures are 2600, 2800 and 3000 K, respectively. The red curves 
are for logg = 4.0 and the blue curves are for logg = 6.0. The spectra were binned to a resolving power of 400 and were normalized separately in each spectral 
band. 



ing from 3.5 to 6.0, and interpolated the measured EW of the 
companion for each respective line. The mean and standard 
deviation of the four lines from this comparison indicate a 
surface gravity of logg = 4.6 ± 0.3, which is relatively low, as 
expec ted and verified for young objects(e.g. Brandeker et al. 
2006). The Kl EWs we measured for the companion are also 
comparable to thos e of other M8-M9 m embers of Upper Sco 
(using spectra from lLodieu et alj200 8) but only a third to half 
of the EWs of M7-M9 fi eld dwarfs of comp arable spectral 
types (using spectra from iMcLean et al. 2003), clearly indi- 
cating low surface gravity and membership in Upper Sco for 
the new companion (see also figure[5]l. 

3.2.2. Atmosphere composition and structure 

Drift-Phoenix provides local gas properties like the gas 
temperature T [K], the gas density p [g cm" 3 ], and the lo- 
cal gas-phase composition, but also dust quantities such as 
the number of dust particles «d [cm -3 ] of mean grain size (a) 
[cm] at each layer of the atmosphere. These quantities are 
needed to evaluate the radiation transfer through the atmo- 
sphere taking into account convective energy transport, hence 
to calculate the synthetic spectrum as shown in Fig. [3] 



The local gas-phase and cloud chemistry is determined by 
the local temperature and pressure which we demonstrate in 
Fig.|6]for all Drift-Phoenix atmosphere models applied in 
Fig. [3] While the low gravity spectra in the 2600-3000 K 
temperature range all provide a reasonable agreement with 
the observed spectrum, figure [6] shows that the correspond- 
ing underlying atmospheric (T,p) profiles are very different. 
This emphasises that medium differences in the spectral en- 
ergy distribution can hide large variation of the atmosphere's 
structure. The models at 3000 K even develop a small temper- 
ature inversion high up in the atmosphere, where H2O is dis- 
sociated, freeing up oxygen and locally enhancing the opacity 
compared to the surrounding layers^ 

Of all the models shown on figure [3] only three do have 
dust in their atmosphere: logg = 4.5 with T e s = 2600 K and 
logg = 6.0 with both T eft = 2600 K and 2800 K. Figure [7] 
shows the mean grain sizes for these models as function of 
local temperature; the lowest temperatures indicate the upper 

8 Such a temperature and pressure inversion does not necesserely invali- 
date the assumption of hydrostatic equilibrium w hich is applied in the DRIFT- 
PHOENIX atmosphere models as was shown in Helling et al. 1 2000) and in 
IAsDlundf!99i3) . 
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FIG. 4.— NIFS spectru m of HIP 78530B (black) compared with the spectrum of USco J155419.99-213543.1 (top), which is an M8 free-floating brown dwarf 
member of Upper Sco l Lodieu et al. 2008), and with a synthetic spectrum from the DRIFT-PHOENIX atmosphere models with 7" c ff = 2800 K, logg = 4.5, and 
solar metallicity (bottom). The spectra were binned to a resolving power of 1000 and were normalized separately in each spectral band. 
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FIG. 5. — Spectrum of HIP 78530B (black) showing the Kl doublets at ~ 1.17 fim and ~ 1 .25 fim. In the top panels, the companion spectrum is compared with 
spectra from the DRIFT-PHOENIX models for T c [f = 2 800 K and logg = 3.5(red), 4.5 (blue), and 5.5 (green). In the bottom panels, the spectrum is compar ed 
with the spectrum of USco J155419.99-213543.1 (red. lLodieu et aU2 008). a young M8 brown dwarf in Upper Sco, and LP412-31 (blue. lMcLean et alj f2003). a 
field M8 dwarf. The shallower Kl lines seen in the companion spectrum relative to the field object indicates lower surface gravity. The spectra were binned to a 
resolving power of ^3000. 
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FIG. 7. — The mean grain sizes across the clouds altitudes in the three dust- 
forming model atmospheres plotted in Fig. [5] The cloud grains remain rather 
small in all models and reside in the optical thin region of the atmosphere. 

The age of the Upper Scorpius association is well con- 
strained at 5 Myr and all stars appear to have formed in 
a burst, over a period of at most 1-2 Myr (Preibis ch et al.l 
120021: ISlesnick et a"i]|2008l) . For an age of 5 ± 1 Myr and the 
above compani on bolometric lumino sity estimate, the evolu- 
tion models of Chabrier et al. (2000) in dicate a mass of ~21- 
26 Mj up , while those of Burrows et al. ( 1997) yield a mass of 
^19-25 Mj up , see figure[8] The ranges of mass quoted reflect 
only the uncertainties on the age and luminosity estimates, but 
of course, larger uncertainty may be present owing to the lack 
of absolute calibration of the evolution models. As visible on 
figure [8] the companion is currently in a phase where its lu- 
minosity is more or less constant with time due to deuterium 
burning; this phase will last until an age of ^15 Myr. Inter- 
estingly, this means that the companion mass estimate is not 
strongly dependent on the age estimate, as is ordinarily the 
case for substellar objects. 

4. DISCUSSION AND CONCLUDING REMARKS 

Due to the very large semi-major axis of the companion 
orbit, its expected orbital motion per year is below the astro- 
metric precision of our data. To estimate the orbital motion, 
we can consider the hypothetical case in which the orbit is cir- 
cular and face-on as seen from Earth, with a semi-major axis 
of 710 AU. Given the system mass of approximately 2.5 M©, 
this yields an orbital period of ^12000 years, which in turn 
gives an angular motion of 0.03° yr" 1 (or a linear motion of 
~2.4 mas yr" 1 ). The angular motion over a decade under 
this assumption, 0.3°, does come out to the same order of 
magnitude as the 2cr differenc e in position ang le between our 
images and those of Kouwenh oven et al.l (120051) . 0.6° ±0.3°. 
However, in reality, the orbital motion is likely to be much 
smaller, as the projected separation is typically smaller than 
the real semi-major axis. In addition, the only way in which 
the angular motion could be larger than 0.3° is if the orbit is 
eccentric and the companion is presently clos e to periastron, 
which is a priori unlikely. As we note in jj3.ll our astrometry 
is internally consistent for the NIRI data but not necessarily 
with respect to other instruments, hence it is plausible that the 
difference with Kouwenhoven et al. is simply spurious. Thus, 
there is no convincing evidence of orbital motion, and no such 
motion is expecte d over the relevant timescales. 

9 |http: //www .ipac.caltech.edu/2mass/releases/allsky/doc/sAqfa4SrfefiWllR pp3ration diagram showing HIP 78530B and 
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FIG. 6. — The temperature-pressure structures for all DRIFT-PHOENIX at- 
mosphere models used in Fig. [5] The T — P profiles vary considerably for 
the three sets of models at different 7/ c ff, although the emerging spectra from 
the three low-gravity atmosphere all provide a reasonable agreement with the 
observed companion spectrum. 



part of the atmosphere. The dust persists at a higher temper- 
ature in the high-gravity models owing to an increased local 
density; however, these models are less relevant for the new 
companion found here, which is young and has low gravity. 
Thus unless the companion is on the low end of our estimated 
T e ff range, it probably does not have dust in its atmosphere. 
However if it is indeed closer to 2600 K and has dust in its at- 
mosphere, the grain sizes would remain rather small, forming 
a haze layer in the upper, optically thin layers of the atmo- 
spheres (see Fig. [TJ- The only objects where haze layers have 
been inferred from tra nsit photometry, h owever, are irradiated 
giant gas planets (see lSing et al.l (120091) ). No direct detection 
of haze on brown dwarfs has been obtained to date. 

3.2.3. Luminosity and mass 

We estimated the bolometric luminosity of the compan- 
ion by using its observed photometry in combination with 
model spectra. We first computed synthetic fluxes in the JHK 
bands using the model spectra, then adjusted those synthetic 
fluxes to the measured values, and finally integrated the scaled 
model spectra over all wavelengths to obtain the total irra- 
diance, which was then converted to bolometric luminosity 
using the primary star distance. The synthetic average flux 
densities in the JHK bands were computed using the rela- 
tive spectral response curves given on the 2MASS project 
webpag^l We repeated t he procedure using b oth the Drift- 
Phoenix and NextGen (Hauschil dt et al.lll999l) model spec- 
tra for ranges of T e g from 2600 K to 3000 K and \ogg from 
3.5 to 5.0. We obtained a value of log(L/L Q ) = -2.55 ±0.13. 
The uncertainty quoted reflects the range of values obtained 
for the different model parameters and accounts for the un- 
certainty on the distance of the primary star. 
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FIG. 8. — Calculated bolometric luminosi ty of the companion ( black p oint) 
compared with ev olutionary tra cks from [Baraffe et al. 1 1998, 20Q2) (red 
curves) and Burrows et al. 1 1997) (blue curves). The curves are labeled with 
their mass expressed in units of solar mass. 

other known brown dwarf and planetary companions to stars 
is shown in Figure [9] With a mass of —23 Mj up and pro- 
jected separation of —700 AU from its massive 2.5 M Q pri- 
mary, HIP 78530B lies at the lower boundary of all known 
companions with separations larger than 100 AU. It is un- 
clear whether these companions, whose mass ratios over- 
lap with those of both more massive brown dwarf compan- 
ions and bona fide planets, formed in a planet-like or in a 
stellar-like manner. Given the size of their orbits, it ap- 
pears unlikely that they formed in situ in a planet-like man- 
ner - either t hrough the collapse of a gravitationally unsta- 
ble disk (e.g. Voro bvqv & Basul 120101) or by the core accre- 
tion mechanism ("e.g. IPollack et al.l 1 1996b . Indeed for both 
scenarios the disks would need to be unusually large and 
massive, but even then, the formation timescale by core ac- 
cretion would be prohibitively long while the efficiency of 
disk instab ility to produce such companions is hig hly un- 
certain dVorobvov & Basul l2010t iKratter et al.ll2010t) . Nev- 
ertheless, it is possible that these companions did form in 
a planet-like manner much closer to the star, but were sub- 
sequ ently kicked outward through gravitation al interactions 
(e.g. IScharf & Menoull2009t IVeras et alJ 120091). I n this sce- 
nario, according to the simulations of IVeras et al.l (120091) . the 
timescale for instabilities to develop and send planets on large 
orbits may be quite short (0.01-1 Myr), with significant dy- 
namical evolution occurring within the first few Myr. How- 
ever, owing to further evolution, planets quickly scattered on 
large orbits are likely to be ejected from the system after a 
few tens of Myr. Alternatively, these wide and low mass ra- 
tio companions could form like stellar binaries, through the 
fragmentation of a pre-stellar core (e.g lBatg|2009t iBate et all 
2003). Numerical simulations indicate that very low mass ra- 
tio companions can indeed be produced by this process, albeit 
rarely, and that they preferentially have large separations and 
are generally found in high order multiple systems. Future ob- 
servations may provide constraints to exclude or support the 
various formation possibilities of these wide low mass com- 
panions. For instance, if they formed in a planet-like man- 
ner closer to the star and were subsequently kicked outward 
through gravitational interactions, then it would be expected 
that additional companions of similar mass or even heavier 
should be present in the systems, at smaller separations. 

In addition, as mentioned above, distant companions pro- 
duced through this process would likely be eventually ejected 
from the system, such that they should be found more fre- 
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FIG. 9. — Mass ratio as a function of separation for various substellar com- 
panions to stars. The red filled circle is HIP 78530B and the red filled di- 
amond is 1RXS J160929.1-210524b. The blue filled circles are all other 
directly imaged low-ma ss substellar companion s (< 25 Mj u „) to stars from 
the compilation given in Lafreniere et al. (2010), augmented with the newly- 
found companion Ross 458C (Goldman et al. 2010; Scholz 2010). The 
large black filled circles are mo re massive directly imaged substellar com- 
panions from the compilation of Zuckerman & Song 12009). The triangles, 
sq uare and down arrow are H R 8799bcde IMarois et at , 2008, 2010), /3 Pic 
b (Lagrange et al. 2009, 2010) and Fomalhaut b iKalas et al. 2008), respec- 
tively. The small black circles indicate planets found by the radial veloc- 
ity and microlensing techniques from the Extrasolar Planets Encyclopaedia 
I http://exoplanet.eu/ 1. 

quently in star-forming regions than around older stars. In 
our overall survey of Upper Scorpius that led to the discovery 
of HIP 78530B and 1RXS J1609-2105b, we observed 91 stars 
(masses 0.15-5 M Q ), so taken at face value, our two discov- 
eries imply that companions with mass ratios below 0.01 and 
separations of hundreds of AUs exist in 2.2^ g% (95% credi- 
bility) of stellar systems. Considering that we did not achieve 
the same sensitivity for all targets as well as our incomplete- 
ness to lower mass ratio companions, this number is only a 
lower limit. The statistics are not yet sufficient to tell whether 
there ex ists a difference for old er objects. For example the 
study of lLafreniereetal.1 (120071) . which targeted —200 Myr- 
old GKM stars, enabled placing upper limits of —6% (95% 
credibility) for c ompanions of —10 M j up (i.e. mass ratio 
-0.01 ); see also IChauvin et al.l d2010l) and iNielsen & Close! 
(2010). Improving the statistics for both young and older 
systems would thus be a good way to investigate further the 
importance of gravitational scattering in accounting for dis- 
tant companions. This approach is valid as long as the in- 
ternal dynamics of the multiple planets dominates over in- 
teractions with other stars or giant molecular clouds as the 
system travels through the galaxy, as these interactions can 
also strip out wide compa nions. This can be roughly checked 
using the results of Weinberg et al.l (119871) . Based on their 
figure 2, with a/M tot — 0.0014 pc M _1 , the disruption life- 
time of HIP 78530AB due to encounters with stars and giant 
molecular clouds is larger than 10 Gyr; this is also the case 
for 1RXS J160929.1-210524Ab. In other words, stellar and 
molecular clouds encounters have little impact on the evolu- 
tion of companions such as the ones we have found. 

The new companion reported in this paper joins a growing 
list of low-mass substellar companions (<25 Mj up ) in wide 
orbit (>100 AU) around stars, wh ich now counts about ten 
objects (see Lafreni ere et al.l (120101) for a recent compilation, 
see also Fig. |9). These companions are found around pri- 
maries covering a wide range of masses and even around pri- 
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maries that are themselves binaries. Thus, while the statistics 
on their frequency may not be very accurate, it seems likely 
that wide low-mass substellar companions are not an unusual 
outcome of the star formation process, yet they remain hard to 
explain within current theoretical frameworks. Over the next 
few years, additional searches for new companions and con- 
tinued efforts to characterize the known ones should allow us 
to make good progress toward explaining their formation. 
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